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C$ ■ The quon algebra gives a description of particles, "quons," that are neither 

fermions nor bosons. The parameter q attached to a quon labels a smooth 
interpolation between bosons, for which q = +1, and fermions, for which 
q = —1. Wigner and Ehrenfest and Oppenheimer showed that a composite 
system of identical bosons and fermions is a fermion if it contains an odd 
number of fermions and is a boson otherwise. Here we generalize this result 
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to composite systems of identical quons. We find q CO mposite = Constituent f° r 
a system of n identical quons. This result reduces to the earlier result for 
bosons and fermions. Using this generalization we find bounds on possible 
violations of the Pauli exclusion principle for nucleons and quarks based on 
such bounds for nuclei. 
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The celebrated spin-statistics theorem JIJ gives the spin-statistics connection 
that lies at the heart of the quantum mechanics of many-body systems: spin one- 
half particles are fermions and integer spin particles are bosons. Wigner || and 
Ehrenfest and Oppenheimer found the fundamental result that composite systems 
of identical bosons and fermions are bosons, unless they have an odd number of 
fermions, in which case they are fermions. Given the usual rules for addition of 
angular momenta, the Wigner-Oppenheimer-Ehrenfest (WEO) result extends the 
spin-statistics connection to composite systems. The WEO result requires that all 
the composites must be in the same internal quantum state and that the interactions 
between the composites must be sufficiently weak not to change the internal states 
of the composites. In this letter we generalize the WEO result to cases where the 
bound states are composed of identical particles that need not be either bosons or 
fermions. 

Recently one of us introduced a formalism, using the quon algebra, that can 
describe small violations of the usual spin-statistics connection [3J. In this formalism 
identical particles can occur in states associated with a mixture of different repre- 
sentations of the symmetric group, with the weights of the mixture dependent on 
a parameter q. To introduce the quon algebra, recall that bosons are described by 
the algebra of commutators of the creation and annihilation operators 

a{k)a\l) -a\l)a(k) =6{k,l), (1) 

and fermions are described by the algebra of anticommutators of the creation and 
annihilation operators, 

a(k)a\l) + a\l)a(k) = 8(k,l). (2) 

The quon algebra is the convex sum of these two expressions, where Eq.(|l]) is mul- 
tiplied by (1 + q)/2, Eq.(|D is multiplied by (1 — q)/2, and the equations are added. 
For the sum to 

be convex q must be between —1 and 1. Several authors have shown that when 
q is in this range all the squared norms of states are positive. || Thus the quon 
algebra is 

a(k)a\l) -qa\l)a{k) =8(k,l). (3) 
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Just as in the case of Bose or Fermi statistics we choose the Fock-like representation 
in which 

a(k)\0) = 0. (4) 

For q = 1 we recover Bose statistics and for q = — 1 we get Fermi statistics. The 
probabilities for two identical particles that are quons to be in the symmetric or the 
antisymmetric state can be found by calculating the following matrix element using 
the rules just given, 

(a*(k 1 )J(k 2 )\0),a%)a%)\0)) = (5) 
5(h, h)5(k 2 , l 2 ) + q5(h, l 2 )8(k 2 , h) = (6) 

^-[5(h, h)S(k 2 , l 2 ) + S(k 1 , l 2 )5(k 2 , h)\ + ^l[S(k 1 , h)S(k 2 , l 2 ) - S(h, l 2 )S(k 2 , h)\, 

(7) 

where we recognize the quantities in the square brackets as the scalar products of 
the symmetric and antisymmetric states of two particles. For systems with more 
than two identical particles there will also be many- dimensional representations of 
the symmetric group. These occur with g-dependent probabilities. For q — > 1 (—1) 
the more symmetric (antisymmetric) representations are more heavily weighted. To 
simplify notation we assume that spin and/or isospin variables and other variables 
are included in the space or momentum variables. 

We want to find a result for the quon statistics of composite states analo- 
gous to that of Wigner and of Ehrenfest and Oppenheimer under the corresponding 
conditions for the internal states and their interactions. First, we note that in gen- 
erating the Fock-like states from the vacuum the order of operators is significant in 
a composite state of quons; in fact for an n-body composite all n! states with per- 
mutations of the orders of the creation operators are linearly independent (provided 
all the operators carry different quantum numbers). This occurs because quons 
can occupy states associated with all the representations of the symmetric group. 
Secondly, we point out that a bound state of quons will be in a single irreducible 
representation, r, of the symmetric group and thus its wave function will have a 
linear combination of terms with the creation operators in permuted orders. Since 
matter is made of electrons and nucleons which are both fermions or, if quons, very 
close to fermions, the most relevant case is when the quons are in an antisymmetric 
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state, r = a, i.e., have q close to —1. Our purpose here is to find the statistical 
behavior of many-particle systems of such bound states. To our surprise, however, 
the statistical behavior of the many-quon bound states is independent of r. The 
only way the representation r enters is in the normalization of the individual bound 
states. To make this clear we first consider a single bound state. 

For a nonrelativistic theory we can represent the bound-state creation operator 
in terms of the creation operators for the constituents as 

/n 
I] d 3 yi f^ (x-y ir --,x-y n )J2c r (P): a\y Pl ) ■ ■ ■ a\y Pn ) : (8) 
i=l P 

or, in momentum space, as 

/n n 
J] d 3 h5(p - J2 ki)f^(k lt ■ ■ ■ , k n ) c r {P) : fitfe) • ■ ■ ~a\k Pn ) : . 
i=i l p 

(9) 

The sums run over the n\ permutations of the symmetric group S n and our permu- 
tations are place permutations rather than label permutations, since place permuta- 
tions are defined even when some of the quantum numbers carried by the operators 
are the same. The coefficients c r (P) pick a sum of products of operators that are in 
the representation r of S n . 

To normalize the bound state we require 

(6 t (a:)|0>,6 t (ic')|0> = S(x-x'). (10) 

To avoid repetition we only discuss the situation in x-space; the calculations in 
momentum space are similar. The normalization condition has the form 

\N r (q)\ 2 V r (q) = 1, (11) 

where V r depends on n (it has degree n(n — l)/2 in q) and also depends on p n >. 
The g-dependence comes from the scalar product 

{■.a\ypi)---a\y Pn ) : |0>, : a\y PV ) ■ ■ ■ a\y Pn >) : |0». (12) 

The graphical rules given in [|J allow calculation of this scalar product. Place the 
operators on one side of the scalar product on one line and those on the other side 
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of the product on a similar line below the first. Each contribution to the vacuum 
matrix element comes from a product of contractions in which a creation operator on 
the top line is paired with a creation operator on the line below to give a 5-function. 
Each set of contractions is associated with a permutation R that takes 

the operators in their original order on the top line and rearranges them to 
the order of the operators with which they are contracted. The power of q associ- 
ated with the product of (^-functions coming from a given set of contractions is the 
inversion number i(R) of this permutation, which equals the minimum number of 
crossings of these lines. Since each operator on one line can pair with every operator 
in the other line, there will be n! terms. This rule is a generalization of Wick's 
theorem in which the usual plus or minus signs are replaced by q l ( R \ 

We now calculate the analog of Eq.(^) with the composite fet's replacing the 
at's. The essential issue in the calculation concerns the algebra of the creation oper- 
ators; the wavefunctions just play the role of deciding which operators are connected 
with a given bound state. Thus only the operator part of the calculation is rele- 
vant. There will be (n!) 2 terms in which the order of the bound states on the two 
sides of the scalar product is preserved but the order of the operators in each of the 
two bound states is permuted in n\ ways. These terms will give two factors of the 
polynomial V r that we found in the normalization of the bound state, one for each 
bound state. The normalization factors obey Eq. flTl]) so they exactly cancel V r . 

The contributions to the scalar product when the bound states are inter- 
changed have a relative factor of q™ 2 in addition to the factor of V 2 just discussed. 
To see this consider the case in which the bound states are interchanged, but the 
order of the operators inside each bound state is preserved. Draw lines from each 
operator in the left bound state on the top line to each operator in the right bound 
state on the bottom line preserving the order of the operators in the bound states. 
Do the same from each operator in the right bound state on the top line to the 
left bound state on the bottom line. The n 2 crossings give the relative factor of q n 
between the terms that are interchanged and those that are not. A similar argument 
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shows that there is a relative factor of q n for all the terms. 

The remaining (2n)! — 2(n\) 2 terms correspond to cases in which creation op- 
erators in one bound state on one side of the scalar product contract with creation 
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operators in both bound states on the other side of the scalar product. These terms 
don't correspond to interchange of the bound state just as they do not in the usual 
Bose and Fermi cases. Here and below we can drop these terms provided we as- 
sume, as in the usual case for the statistics of bound states, that the interaction 
between the bound states is negligible compared to the interaction that binds the 
constituents in each bound state. Thus we have established that the analog of Eq.(^) 
for the bound states is 

(6 t (pi)6 t (P2)|0>,6+(ri)6+(r 2 )|0))= (13) 

6(pi,r 1 )6(p 2 ,r 2 ) + q n 5(pi, r 2 )6(p 2 , n) = (14) 
]_ + q n 2 ]_ _ q n 2 

— 2 — [S{pi,r 1 )5{p2,r2)+5{pi,r 2 )5{p2,r 1 )]-\ — [<J(pi, r 1 )6(p 2 , r 2 )-6{p 1 , r 2 )5(p 2) n)\. 

(15) 

We can now represent the n\ lines that connect a bound state on the left hand 
side of a scalar product to a bound state on the right hand side by a superline. When 
we use our generalized Wick's theorem to calculate matrix elements of the bound 
state operators we find the same result as for a single quon operator except that the 
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number of crossings of superlines replaces the number of crossings of lines so that q n 
replaces q. Thus bound states of n quons with parameter q have quon statistics with 
parameter q n2 . || |7j This generalization of the Wigner-Ehrenfest-Oppenheimer rule 
is the main result of this letter. 

Since Bose and Fermi statistics correspond to q = 1 and q = — 1, respectively, 
we expect the Wigner-Ehrenfest-Oppenheimer rule to emerge in the limit q — > ±1. 
Indeed this is true since n and n 2 are even or odd together. 

The composite statistics rule is important in establishing the statistical behav- 
ior of identical nuclei in molecules for example. Similar issues arise in the collective 
behavior of Cooper pairs in BCS-type superconductivity [T5j and superfluid behav- 
ior in 3 He [ 16fl . The crucial condition for this rule, as pointed out by Ehrenfest and 



Oppenheimer, is 

that all the composites must be in the same internal state and the interactions 
between the composites must be sufficiently weak so as not to change the internal 
states. In other words the interaction energies among the composites must be small 

6 



compared to the energy separation of the internal states. This condition is obviously 
well satisfied for nuclei in molecules where the interactions are on the order of 
electron volts and the internal energy differences are on the order of MeV. On the 
other hand nitrogen atoms taken individually are fermions (seven electrons, seven 
protons and seven neutrons for 14 N) but the N 2 molecule is permutation symmetric 
under the interchange of the two 1A N nuclei; the statistics of the individual atoms 
is no 

longer relevant once the atoms interact to form the molecule. 

We now apply the composite statistics rule to nuclei. We don't expect exper- 
iments using bound states of many electrons, i.e., atoms, to improve the extremely 
high-precision bound for electrons due to E. Ramberg and G.A. Snow ||. On the 
other hand experiments with nuclei may give significant bounds on the statistics of 
nucleons, for which there are at present no direct precision experimental tests. For 
nuclei that are bound states of A almost-Fermi nucleons the result is immediate. If 
A is even so the nucleus is close to Bose with parameter 1 > g# = 1 — e, then the 
constituent nucleons are quons with parameter — 1 < q nU cieon = —I + e/A 2 . If A is 
odd so the nucleus is close to Fermi with parameter — 1 < qp = — 1 + e, then the 
constituent nucleons are quons with parameter — 1 < q n udeon = —1 + A 2 . Thus 
the bound on the deviation of the nucleons from Fermi statistics is improved over 
the bound on the deviation of the nucleus by a factor of 1/A 2 . This analysis can be 
extended to quarks bound in nucleons or other baryons, for which q\ uark = qbaryon- 

Empirical bounds on violations of statistics found up to 1989 are reviewed in 
Pj. More recent experiments have set stringent upper limits on the violation of the 



spin-statistics connection for electrons [H, [L0| and 16 nuclei O, |12| . DeMille and 



Derr have found a bound on violation of Bose statistics for photons ||13|| . From a 
sensitive C w O<? 



molecular spectroscopy experiment Modugno, Ingusicio, and Tino |L4J have 
been able to show that the probability of finding the two 16 nuclei (spin 0) in 
carbon dioxide in a permutation antisymmetric state is less than 5 x 10 -9 . 

The 16 nucleus is of course a composite system of eight protons and eight 
neutrons. Applying the results derived above, we can interpret the experimental 
limit on spin-statistics violations for 16 as setting a limit of 1 + q nU cieon < 2 x 10~ n 
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for nucleons. We believe this is the first precision limit on the violation of the spin- 
statistics connection for nucleons linked directly to an experimental measurement. 
Other limits on violations of the Exclusion Principle for nucleons have relied on 
arguments based on models of the solar p-p cycle |I7| (from which a limit of about 
10 -15 is claimed) or arguments based on models for isotopic elemental production 
in supernovae core collapses fll8 |. In the latter case there are many uncertainties in 



the model for the collapse and the precision mass spectrometry searches for unusual 
isotopes have yet to be 

done. Extending the analysis to quarks within the model of nucleons as bound 
states of 3 quarks, yields 1 + q quar k < 2 x 1CT 12 for quarks. 

There is yet another way to extend existing experimental limits. In a recent 



paper Greenberg and Hilborn [19| have argued that an experimental limit on the 
spin-statistics violations for one type of particle, say electrons, can be used to set 
a limit on spin-statistics violations for particles with which the original particles 
interact. For example it is possible to extend the stringent limit 1 + q < 10 -26 for a 
spin-statistics violation for 

electrons from the Ramberg-Snow experiment [§] to set a limit on a possible 
spin-statistics connection violation for photons. 

The theorem proved in this paper provides the quon generalization of the 

standard rule for the statistics of composite systems (fermion behavior for an 
odd number of fermions; otherwise boson behavior occurs). In accord with the 
original work by Wigner and by Ehrenfest and Oppenheimer, we have shown that 

2 

the generalized rule q compos ite = qconstituent holds as long as the composites are all in 
the same internal state and the interaction energies among the composites can be 
ignored compared to the differences in internal state energies (i.e. the composites 
are tightly bound). 
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